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Abstract. We show that the spin-charge separation predicted for correlated fermions 
in one dimension, could be observed using polarized photons propagating in a nonlinear 
optical waveguide. Using coherent control techniques and employing a cold atom 
ensemble interacting with the photons, large nonlinearities in the single photon level 
can be achieved. We show that the latter can allow for the simulation of a strongly 
interacting gas, which is made of stationary dark-state polaritons of two species and 
then shown to form a Luttinger liquid of effective fermions for the right regime 
of interactions. The system can be tuned optically to the relevant regime where 
the spin-charge separation is expected to occur. The characteristic features of the 
separation as demonstrated in the different spin and charge densities and velocities can 
be efficiently detected via optical measurements of the emitted photons with current 
optical technologies. 
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1. Introduction 

1.1. Spin-charge separation and quantum simulators 

One-dimensional (ID) physical systems have attracted much attention due to their novel 
and sometimes spectacular features. Unlike two- or three-dimensional systems, the 
physics of ID Lieb-Liniger model pj is well captured by the Luttinger liquid theory in 
the low-energy domain |2J. In ID Luttinger liquid theory, collective excitations rather 
than single excitations appear due to the tight transverse confinement forcing particles 
to move along one direction and thus converting any individual motion into a collective 
one. The collective modes of spin and charge in the ID electronic gases surprisingly can 
be shown to propagate with different velocities, i. e., the spin and the charge separate into 
spinons and holons [2J [3, HI E] • ^ n experiments, the observation of spin-charge separation 
is however challenging - the control of interactions is still elusive and no distinct features 
of separation are obtained, although several seminal works have been made using copper- 
oxide chain compounds SrCu02 |SJ, metallic chains [7J, superconductors [8], and more 
recently, GaAs/AlGaAs heterostructures [9] [TO] ITT]. 

At the same time, works on artificially engineered quantum optical systems in 
which many-body effects could be reproduced in well controlled environments have 
recently emerged. To observe spin-charge separation in cold atoms, experimental 
proposals involving bosonic and fermionic species have been suggested [121 03] EH EE] . 
However, the challenges in the trapping and cooling of fermionic gases, especially 
the lack of necessary individual accessibility and measurement of correlations in 
general, make current results inconclusive. Strongly correlated photons and polaritons, 
as hybrid light-matter quantum simulators, have been recently proposed. Initially 
using coupled resonator implementations [IB] . Mott transitions [XT] [JSJ [19] and then 
effective spin models and Fractional Hall states of light were shown to be possible 
[201 I2T1 [221 [231 EMI [25] . More recently, by employing hollow-core optical waveguides 
filled with cold atom ensembles and using slow-light techiques [26], [27], [28] , it was shown 
that it is possible to prepare a Tonks gas of photons [29]. Using two atomic species, a 
two-component Lieb-Liniger model has also been recently suggested [30] . Compared to 
cold-atom proposals, photonic proposals should allow for more direct measurements of 
local observables and correlation functions of the emitted photon states. 

In this article, we consider a novel scheme involving two oppositely circularly 
polarized quantum beams and single atomic species to simulate a two-component 
interacting gas in ID. We note the difference to an earlier scheme, where two quantum 
fields with different frequencies interacting with two species of four-level atoms were 
employed [30]. The current proposal has the distinct advantage over in two major 
ways. Firstly in terms of more efficient detection of the correlation in output polarized 
states and secondly in terms of the loading and preparation process into the fiber which 
requires handling a single atomic species rather than two. 

The paper is organized as follows. We first review the basics of the single- and two- 
component bosonic Lieb-Liniger model [1] , and especially its mapping to Luttinger liquid 
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theory. We then describe in detail the preparation of a polaritonic Lieb-Liniger model 
with two bosonic components in the waveguide employing slow-light techniques. Finally 
we discuss the necessary range of optical parameters in order to drive the system to the 
relevant spin-charge separation regime. We conclude by analyzing how the effective 
photonic spin and charge densities and velocities can be probed. This is done by 
releasing the trapped polaritons into outgoing photons where optical measurements can 
reconstruct the characteristic functions of the effect. 

1.2. From Lieb-Liniger bosons to Luttinger liquids: The basics of Spin- Charge 
Separation revisited 

One of the most famous ID models is the Lieb-Liniger model [I], which describes iV 
bosons with a Dirac-delta interaction as 



where m is the mass and g is the interaction strength. Although the Lieb-Liniger 
model is exactly solvable through the Bethe ansatz approach, it is still generally very 
difficult to extract correlation functions from the solutions at any interaction regime. 
Luttinger liquid theory on the other hand can give the low energy universal description 
of the Lieb-Liniger model at low temperatures. In the Luttinger liquid phase, the low 
energy excitations are no longer single but collective modes with linear dispersion. The 
confinement of interactions between particles in ID forces any individual motion to affect 
the collective system. The description of the dynamics in terms of collective bosonic 
fields is called 'bosonization' approach and we briefly present it in the following. 

Assume N particles described by ^(r) = ip( z )¥o(. r ±)> where the particles move along 
one direction, say z direction. Strong confinement is applied to the transverse direction 
r ± = {x, y}, allowing only for the lowest energy transverse quantum state (po( r ±) to be 
considered. The Lieb-Liniger model of the parallel component reads [21 EJ, SI EJ 



with p (z) = ip^ip the particle density and 9 (z) the phase operator. As described in [2], 
the phase and density fields are canonically conjugated fields, 




(1) 




(2) 




(3) 




(4) 



and the collective density operator can be expressed as 




(5) 



m=— oo 
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If p(z) varies slowly with z, we can retain the lowest frequency component for m = 
and write 

pO) - Po + -d z <f){z), (6) 

7T 

where the fields 9(z) and ^d z (f>(z') are canonically conjugated. 

The Lieb-Liniger Hamiltonian (2) for the lowest component ip(z) ~ e~ ie<yZ ^ p 1 ^ 2 can 
be mapped to a Luttinger liquid with Hamiltonian (see [21 El IH 

RS = J S [ ^ Sl( ^ )2 + ]^ ( ^ )2] ' (7) 
where all the interaction effects are encoded into two effective parameters: the 
propagation velocity of density disturbances v and the so-called Luttinger parameter 
K sl controlling the long-distance decay of correlations. 

Going beyond the simple case with single- component bosonic system, interesting 
physics can be obtained by mixing two bosonic components or by considering two 
internal degrees of freedom of bosons, which is analogous to assigning a "spin" 1/2 to 
bosons. The two-component Lieb-Liniger model in this case would read: 

H*= fdzJ2 \^-d z iPt(z)d z ip s (z) + y /£(*)] + / dz X nP^) Pi (z), (8) 

where m s is the mass and Xs and xu are the intra- and interspecies interaction with 
s =t, 4- representing the two spins. Following the literature [21 El IH E], the charge- and 
spin- related fields can be defined as: 

^ charge ^j-^ i "spin ' ^charge ' 'rspin ' V / 

The Hamiltonian separates into two parts (more details in [2]) and reads as: = 
-^charge + H spin . The charge part is given by 

-^charge / 7) [^charge -^charge {&x ^charge) ~\~ ~T7 (^x0charge) ] (l 0) 

J -^charge 

and the spin part is defined as 



/dz 2 ^spin 2 

~Z [Wspin-^spin(^2:^spin) ~T7 (^:r0spin) ] 
Zn A S pi n 

+ 2 XUPo COS (v^spin) (11) 

under the separation conditions 

Xt = Xi, = • 12 

Here po = Po,t + Pcu and po,s are the densities for two species. With x — Xs, 
u = u s = a/ Po,sXs/ m s> and K = K s = y7r 2 p ,s/(^sXs)j the charge and spin velocities 
and Luttinger parameters are u c har g e, S pm = Uy/l ± Xn/x, -^char g e, sp in = K/ >/T± Xn/X- 
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2. Photonic Spin-Charge Separation in Nonlinear Optical Waveguides 

2.1. The Optical Waveguide Setup 

Our proposal is based on exploiting the available huge photonic nonlinearities possible 
to generate in specific quantum optical setups. More specifically, we envisage the use 
of a highly nonlinear waveguide where the necessary nonlinearity will emerge through 
the strong interaction of the propagating photons to existing emitters in the waveguide. 
Recent experiments have developed two similar setups in this direction, both capable 
of implementing our proposal with either current or near future platforms. In these 
experiments, cold atomic ensembles are brought close to the surface of a tapered fiber 
[3TJ [32] or are loaded inside the core of a hollow-core waveguide [33j |3U |35j [36] as shown 
in figure [I] (a). The available optical nonlinearity based on the Electromagnetically 
Induced Transparency (EIT) effect can be used as we will show to create situations 
where the trapped photons obey Lieb-Liniger physics. 

The process to generate the strongly correlated states of photons is as follows: First, 
laser-cooled atoms exhibiting a multiple atomic-level structure shown in figure [l] (b) are 
moved into position so they will interact strongly with incident quantum light fields. 
Initially, resonant to the corresponding transitions, two optical pulses with opposite 
polarizations E^ }+ (z,t) and E^ + (z,t) are sent in from one direction, say the left side. 
They are injected into the waveguide with the co-propagating classical control fields 
and Q± i+ (t) initially turned on. As soon as the two quantum pulses completely 
enter into the waveguide, the classical fields Q St+ are adiabatically turned off, converting 
E s , + into coherent atomic excitations as in usual slow-light experiments for s =t, \,. 
We then adiabatically switch on both f2 s>+ and f2 Sj _ from two sides. The probe pulses 
become trapped due to the effective Bragg scattering from the stationary classical waves 
as analyzed in [2EJ E7J EE] • At this stage the pulses are noninteracting with the photons 
expanding freely due to the dispersion. By slowly shifting the (i-levels, the effective 
masses can be kept constant whereas the effective intra- and interspecies repulsions are 
increased. This drives the system into a strongly interacting regime. This dynamic 
evolution is possible by keeping for example the corresponding photon detunings A s 
constant while shifting the d-level. Once this correlated state is achieved, the fields - 
for example Q s .+ - from the pair of control fields that trap polaritons, are slowly turned 
off. This will release the corresponding quasi-particles by turning them to propagating 
photons which will then exit the fiber. As all correlations established in the previous 
step are retained, these wave packets comprise of two separated effective charge and 
spin density waves. 

2.2. Realizing a two-component Lieb-Liniger model of polarized photons 
The system described above and shown in figure [T] obeys the Hamiltonian: 
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+ V / 27r(^ s CT 6)S . a + 9ss'^d,s,s';c,s') 
s' 

x (E Sj+ e ih ^ sZ + E Si _e~ ik ^ z ^ 

+ {Q s , + e ikc - z + Q s ^e~ ikc - Z ) a CjS;M + h.c.}, (13) 

where s, s' =4, 4 and \d, f, 4) = |<i, 4, T)- The continuous collective atomic spin operators 
a mv = o^ v {z,t) describe the averages of the flip operators \fj) (y\ over atoms in a small 
region around z. The density of atoms is n s z and g s , g ss i are the coupling strengths 
between the quantum fields and atoms, while A s and A ss i are one-photon detunings 
from the corresponding transitions. For simplicity, we assume that g s = g ss > = g. 
Furthermore, we label the two quantum and two classical fields with frequencies uq >s 
and wc, s and wave vectors &q iS and kc, s , respectively. Both quantum fields E^ }+ (z,t) 
and E± t+ (z,t) drive four possible atomic transitions. The fields E Sj ±(z,t) are detuned 
by A s from the transition \a) — > \b,s) and by A ss from \c,s) — > \d,s,s). E St ±(z,t) also 
drive the transitions from |c, s) — > \d, s, s) with detuning A sS . Here s =4,4 and s^s. 
Finally, the applied classical control beams with Rabi frequencies n s ±(t) couple to both 
atoms and drive the transitions \b,s) — > \c,s). 

The evolution of the slowly-varying quantum operators E s> ± are given by four 
Maxwell-Bloch (MB) equations 

(d t + vd z ) E Sj± = iV27in s z g (cr a;fejSi± + o- c ,s;d,s,s,± + (T c ,s;d,s,s,±) (14) 
with four levels of the s-th atoms denoted as \a, s), \b,s) , |c, s) and \d,s,s). When 



writing down the MB equations (14), we have introduced the slowly- varying collective 
operators 

<r K u = <w(z,t)e ife ^* + a^(z,t)e~ ik ^ z , (15) 

and v is the velocity of quantum fields in an empty waveguide [§} In the derivation of 
equations of motion we assume the Rabi frequencies of the control fields to be slowly 
varied. The slow-light polariton operators are defined as 

ty St ± = cos 9 s E Si ± - sin 9 s ^2nn s z o- C)S . a , (16) 

where tan# s = gy/2irn z /Q s . For stationary polaritons we have assumed that the 
amplitudes of the counterpropagating classical fields are equal, i.e., Q Sj ± = Q s . In 
the limit when the excitations are mostly in spin-wave form, i.e., sm8 s ~ 1, and since 
o~ c ,s;a = —gE St ±/Q s , the polariton operators become 

* Si± = v/2^^4 ± . (17) 

Setting ^ s = (^ Sj+ + * s -)/2 and A s = (fy s>+ - * s ,-)/2 as the symmetric and 
antisymmetric combinations of the two polaritons, we derive the equations of motion 

§ We would like to highlight here the relative simplicity of the above evolution equation compared 
to the one we considered in |30j where extra phase terms have to be involved due to the existence of 
two-atomic-species different frequencies on the incident quantum fields. 
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Figure 1. The model setup under consideration. In a nonlinear fiber ( a hollow core 
version is shown here but a tapered fiber approach could also be used) , the cold atoms 
are interacting with two quantum light fields E s (red and blue arrowlines) and two pairs 
of classical fields fi fll ± (yellow and green arrowlines), where ± denotes the forward or 
backward propagation direction. The studied atomic level diagram and possible atomic 
transitions driven by two oppositely circularly polarized quantum pulses E s and two 
control beams il s is shown in (b) . Appropriate tuning the couplings of light fields to 
the corresponding atomic transitions, forces the trapped polaritons to behave as an 
effective ID quantum Luttinger liquid and reach the spin-charge separation regime. 
Coherently transferring the polaritons' correlations to propagating light pulses and 
allowing them to exit the fiber, provide for the efficient measurement of both the 
dynamics of the propagation of the effective spin and charge quasiparticles and the 
effective spin and charge velocities characteristic of the effect taking place. 



for the polariton combinations A s : 

A 



d^ 8 + vd z A s = -7T tan 2 9 s d t ^s ~ (2***, + A\A.) * 4 



- h¥^s + At As) ^ s + noise, 

d t A s + vd z V 8 = -i^-n' s A t - + noise. (18) 



The noise terms in Eq. (18) account for the dissipative processes that take place 
during the evolution. Fortunately, for the dark-state polaritons under consideration, as 
long as the spontaneous emission rates T from the states \c,s) and \d, s,s') are much 
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less than the detunings |A sa /|, the losses in the timescales of interest are not significant 
and thus can be neglected as discussed in [26], 123 1211 I2H1 [30]. Assuming optical depth 
of a few thousand and a large ratio between the density of atoms to the density of 
photons n s z /n s vh ~ 10 4 , the antisymmetric combinations A|- and A± can be adiabatically 
eliminated from the equations of motion for the polaritons and moreover, the nonlinear 



terms like ^>\^ S A S and Ajy4 s \l/ S are negligible. In this regime, Eq. (18) simplifies to a 



nonlinear Schrodinger Eq. (19) for polaritons which reads: 

r 



2A s za 



ID 



tr: 



A,, 



+ 



A, B s 



(19) 



which is related to an effective two-component Lieb-Liniger model of polaritons 

1 



H 



dzj2 



—dM(^s(z) + ^pl{z) 



+ / dz Xn Pt(z)Pi(z).(20) 



-Y s 1Tl n s z / (AA s v s ) is the effective mass for s-th polaritons with T 



s 

ID 



Here m s = 

47rg 2 /V the spontaneous emission rate of a single atom into the waveguide modes and 
v s = vVL 2 s / (Trg 2 n s z ) the group velocity of the propagating polaritons. The intraspecies 
repulsion is characterized by Xs = Tid^/Ass and the interspecies repulsion by xu = 

To reach the spin-charge separation regime, we employ the mapping to a Luttinger 
liquid model described in the previous section. For this to be possible, we should first 
of all check the tunability of the relevant parameters to the repulsive regime. This in 
our case (see Eq. (JT9[)) implies tuning A S A SS < which leads to m s x s > 0. Similarly 
A-fA^ > forces the effective masses m^m^ > and A^A^ > tunes x^Xi > 0- These 
conditions are satisfied by tuning the lasers such that: A^ and A^ are negative (positive) 
while at the same time A^ and are positive (negative). 

Apart from the repulsive interaction regime, the separation condition Eq 
needs to be satisfied as well, which in our case means setting: 

Po,t Pcu 



(12) 



(21) 



where the polariton density p ,s equals to the photon density n^ h . The effective charge 
and spin densities are the sum and difference of the two-species polaritonic densities, 
which read as 

Pcharge = Pt + Ph Pspin = P\ ~ Pi (22) 

with p s = (p ,s - ^V0 S ) XL ex P l i2m ( n Po,s z - 0s)]- Keeping only the lowest 
components with m = 0, ±1, the charge and spin density operators in the bosonic 
language can be represented as 



V% r r 

Pcharge = Po ^charge + 2 />0 COS 2k F Z - V20 char 



COS 



^20; 



Pspin 



y/2 







7T 



zS^spin 



+ 2p sin 



2k F z - A/20 charge sin \/20, 



y spm- 



spin; 

(23) 
(24) 
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Figure 2. The Fourier transform of density-density correlations for our polaritonic 
system. It exhibits the characteristic splitting corresponding to the two different 
propagation velocities for the photonic spinons and holons, as calculated in [37 . The 
effective spin and charge velocities are u c hargc = 1 and u sp in = 0.5 (normalized by 

2y/|u), and the Luttinger parameters are ^charge = 0.55 and K sp i n = 1.1. They 
translate in our case at optical depths of OD = 2000, 10 photons in each pulse initially, 
and single-atom co-operativity for each atomic species of 0.4. In (b) a cut of the 3D plot 
at quasimomentum q = 2 is plotted to show the distinct two peaks corresponding to the 
spin and charge velocities. The cross density-density correlations are experimentally 
reconstructed via typical optical measurements on the correlated photon states as they 
exit the fiber. 



Here the first term in p c h ar g e is the average density p = Po,t + Po,i- in our two-species 
photonic system, we set p ,t = Au = \po for each polarization component. The second 
gradient term in p c h a rge an d p sp i n are the density oscillations with zero momentum. The 
third term in p c har g e and second term in p spin are the density fluctuations of the 2fcp 
components [37]. We label 7 S as the ratio of the interaction to the kinetic energies for 
each polariton species 7 S = m s x s /po,s- Combining the two separation conditions in Eq. 
(21) together, one gets ^ = 74.. For x = Xs an d 7 = j s , the velocities and Luttinger 
parameters can be expressed as u = x/s/l an d K = Tr/y^y. As also demonstrated for a 
similar system albeit with one quantum field [29J, 7 here can also be tuned from zero to 
finite to extremely large, corresponding to non-, weak- and strong-correlated regimes, 
which implies a wide tunable range for u and K. 



2.3. Probing of the photonic spinons and holons 

The typical detection of spin-charge separation can occur through dynamically probing 
the time evolution of a single excitation as in cold-atom proposals [121 [131 EH 03] > or 
by measuring the corresponding single-particle spectral function as in condensed matter 
experiments [2],[3llH[5l[6l[71|8l[9l H21 E]- m ° ur case, we propose to extract the charge 
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and spin velocities by measuring the Fourier transform of density-density correlations 
D(uj, q) for energy uj and momentum q. As derived for the two-component system in 
for the 2k-p component of density operator p 2 fc F = 2po cos (2kpz — v^Acharge) cos v^20: 



spin 



(the last term in Eq. (23 )), the Fourier transform D(u,q) of the density-density operator 



(T T p2k F (z, r) p 2kF (0, 0)) is given by 

4vrp2 (D^W+i^in r (1 _ Kchargc/2 _ K spin/2 ) 



D(u,q)= - 



x exp 



-111 



r (-^-charge/ 2 "I" -K"spin/ 2 ) 
1 1 

(-^charge + -Kspin \ _ / 2 2 2\ 
^ 1 J 8 (w - M spin g ) 



v I. ,2 „,2 2| i ^charge/2+-fCspin/2-l l--^ct 
X | W "spiny | "spin 



j-! / -^charge -^charge -^-spin 1 -, -^-charge -^spin 
X o ) o > l 



2 ' 2 2 

-^charge H~ -^-spin -. ^charge ^ ^charge 1 ? \ /nr\ 

2 "spin ^ - <in<? 2 

where T is the gamma function, G is the step function, F\ is the Appell's hypergeometric 
function, and a is a short-distance cutoff. D(u>,q) depends on the velocity oj/q and 
should exhibit two peaks centered around Mcharge? and u spin q [21 El IU IS] • In our photonic 
system, probing of the spinon and holon branches can be done by measuring the 
correlation functions of densities of the fields as they exit, for a specific quasi-momenta 
q. For a clear distinction between the two effective spin and charge peaks, we should set 
our optical detectors around q = 2ti/z , i.e., z apart (z here corresponds to roughly 
the length of the fiber). To give an illustration of the expected behavior, D(u,q) in 
the unit of p\a is plotted in figure [2] with the intra- and interspecies repulsion ratio 
Xt-i/x = 0-6, which in turn tunes the charge and spin velocities to 

2 



^charge = 2w spin = U\ 1 + Xu/X = 2 A/ ~U, (26) 



and Luttinger parameters -Kdiarge = lK spin = \\J\K. We choose w c harge = 1, ^s P in = 0.5 
and -R'charge = 0.55, ifgpin = 11 via tuning u and K wihch require 7 ~ 20. This as 

A- I 

shown in [29] is achieved at optical depths OD = 2000 and roughly N^t = 10 photons 



ph 

initially in each pulse and single-atom co-operativity of 77 = 0.4 [jj] These values are 
for the moment out of the current experimental range where optical depths of a few 
hundred have been achieved, but should not be out of the question in the near to mid- 
term future [211 E21 EE]- In calculating the optical interaction parameters appearing in 
the Hamiltonian Eq. (20), we haven take into account both the linear and nonlinear 
loss mechanisms as layed out in 



Co-operativity here is the ratio of spontaneous emission into the waveguide to total spontaneous 



emission. 
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3. Conclusion 

We have described in detail a strongly correlated photonic scheme to simulate a 
purely fermionic effect, spin-charge separation. In more detail, we have shown that 
polarized photons interacting with a cold atomic ensemble can be made to obey two- 
component Lieb-Liniger physics and even behave as a quantum Luttinger liquid. The 
relevant interactions exhibit the necessary tunability for steering the photons to the 
effective spin-charge separation regime. Efficient observations of the characteristic 
features of the separation using standard quantum optical methods should be feasible 
based on correlations measurements of the outgoing photons which here carry opposite 
polarizations. The current proposal is different from a similar scheme proposed earlier 
by some of us, where two species of atoms were coupled to two quantum fields of two 
different frequencies but of same polarization [30]. Here a single species of atoms is 
shown to suffice in order to induce the required intra- and interspecies interactions, 
which combined with the easier detection of the polarized output states makes this 
approach more feasible. 
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